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^ Effect of a lattice upon an interacting system of electrons in 
1^ two dimensions: Breakdown of scaling and decay of persistent 
cn currents 
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rj-^ When one considers interacting electrons free to move in energy spectrum. What is the carrier density above which 

^ ; an immobile background of positive ions, one can repre- the scale s becomes relevant and ceases to be the scal- 

O , sent the ions by a uniform positive jellium if the electron j^g parameter for the lattice ground state? The answer 

O ■ density is sufHciently small. This uniform jellium gives depends on s and on the two parameters controlling the 

^ [ simply rise to a constant term in the Hamiltonian. One effective Bohr radius 
• . gets a continuum model characterized by two scales: the 

' inter-electron spacing a and the effective Bohr radius a|j. Srh^ 

^ i Simple scaling laws are obtained if one uses the dimen- = ^*q2 ■ (■'■) 
■ - - ' sionless ratio Vs — a/a*^. This continuum approximation 

neglects the discrete character of the lattice of positive the effective mass m* of the carriers and the dielectric 

ions. constant Sr of the medium. 

If one wants to keep the periodic character of the ionic If a*g is large compared to s, the lattice effects are only 

lattice, one has to include a periodic potential instead of a important for small values of r^, for which one has a quan- 

uniform jellium or to use the tight-binding approximation, tum fluid with a deformed Fermi surface. This is a highly 

One obtains a lattice model, where the kinetic energy can quantum weak coupling limit of large carrier densities. If 

be simplified if the hopping terms are restricted to nearest the issue is to study charge crystallization in such sys- 

neighbor ions. A lattice introduces a third scale: the lattice terns, the densities of interest are much lower than those 

spacing s. If s is irrelevant, the lattice model keeps the required to deform the Fermi surface, and the physics can 

same universal scaling than the continuum limit, if one be described in the continuum approximation. If one uses 

uses the combination of lattice parameters which becomes a. L x L lattice model with N particles to study electron 

Ts ~ in the continuum. crystallization for systems of small s/a*^, one has to take 
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Abstract. The ground state of an electron gas is characterized by the interparticle spacing to the effective 
Bohr radius ratio — a/a*g. For polarized electrons on a two dimensional square lattice with Coulomb 
repulsion, we study the threshold value rl below which the lattice spacing s becomes a relevant scale and 
Ta ceases to be the scaling parameter. For systems of small ratios s/a*g, s becomes only relevant at small Vs 
(large densities) where one has a quantum fluid with a deformed Fermi surface. For systems of large s/a*g, s 
plays also a role at large Vs (small densities) where one has a Wigner solid, the lattice limiting its harmonic 
vibrations. The thermodynamic limit of physical systems of different a'g is qualitatively discussed, before 
quantitatively studying the lattice effects occurring at large r^. Using a few particle system, we compare 
exact numerical results obtained with a lattice and analytical perturbative expansions obtained in the 
continuum limit. Three criteria giving similar values for the lattice threshold r* are proposed. The flrst 
one is a delocalization criterion in the Fock basis of lattice site orbitals. The second one uses the persistent 
current which can depend on the interaction in a lattice, while it becomes independent of the interaction 
in the continuum limit. The third one takes into account the limit imposed by the lattice to the harmonic 
vibrations of the electron solid. 

PACS. 71.10.-w Theories and models for many-electron systems - 71.10.Fd Lattice fermion models - 
73.20.Qt Electron sohds 

1 Introduction We study when the low energy spectrum of the lattice 

model can be described by a continuum approximation, 
the lattice effects remaining only important in the high 
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lattice fillings N / 1? and hence ratios above the lattice 
threshold r* , where the lattice model can be described by 
the continuum limit. 

If a*Q is small compared to s, the lattice effects be- 
come important when is large also. This is a limit where 
the lattice physics takes place also in a weakly quantum 
strong coupling regime, at low densities. If one continues 
to increase the density in those systems, one can eventu- 
ally reach the limit usually described by a Hubbard model 
near half filling, where the lattice can give rise to a Mott 
insulator if the interaction is large enough. In contrast to 
the case where is large compared to s, the continuum 
approximation cannot be assumed for studying electron 
crystallization. One has the obvious problem of commen- 
surability between the electron lattice characterizing the 
continuum limit and the ionic lattice. Even if there is com- 
mensurability, there is a remaining limit imposed by the 
lattice to the harmonic vibrations of the Wigner solid. 

Eventually, let us note that for a lattice model, it is im- 
portant to know the electron density below which its low 
energy spectrum begins to exhibit the continuum behavior 
and its universal scaling laws, if one uses the combination 
of lattice parameters which becomes = a/a|j in the 
continuum. 



2 Two dimensional continuum model 

The Hamiltonian ifc describing N polarized electrons of 
mass m free to move on a continuum space of dimension 
d and dielectric constant = 1 contains one body kinetic 
terms, two body interaction terms plus the constant term 
due to the presence of the uniform background of positive 
ions necessary to have charge neutrality. 



.2 W 



Br 



2m 



1 



const, 



i=l 



l<i<j<N 



Measuring the energies in rydbergs {IRy — me*/2?i^) and 
the lengths in units of the radius a of a sphere (circle in 
2d) which encloses on the average one electron, e and m 
being the electronic charge and mass. He becomes 
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-t- const, (3) 



be easily extended to arbitrary dimensions and could in- 
clude the spin degrees of freedom, we restrict the study to 
polarized electrons (spinless fermions) and to the dimen- 
sion d = 2. The ground state (GS) energy in rydbergs per 
particle is given in the Hartree-Fock approximation [J as: 



i?o = ^ + ^ + 0(lnr. 



< 1 



(4) 



with coefficients ho = 2 for the kinetic energy and hi = 
— 1.6972 for the exchange energy. 

In the strong coupling limit (large r^), the volume per 
electron a'' is large compared to a^, and the electrons crys- 
tallize on an hexagonal lattice with weak quantum effects 
(Wigner crystal). As Wigner's original approximation^] 
suggests, the GS energy per particle in rydbergs can be 
expanded in powers of 



1/2 



Eo 



^ + 4 + § + o(.f/^), 



rs » 1 



(5) 



l<i<j<N 



The leading term (oc r~^) is of classical nature (Coulomb 
energy of the lattice of electrons in a continuum back- 
ground of positive charge) while the first correction (cx 

rs ^^'^) is quantum (zero point energy of the harmonic os- 
cillations of the electrons about their lattice points). One 
gets piim^ /o = -2.2122 and /i = 1.628 in two dimen- 
sions. 

Using quantum Monte Carlo methods, the two dimen- 
sional crossover between these two limits has been studied. 
A variational approach and a Green function Monte 
Carlo approach 6 have given a critical ratio « 37 for a 
possible transition separating the quantum fluid from the 
Wigner solid in the continuum. However, the well-known 
sign problem of the Monte Carlo methods requires to im- 
pose the nodal structures of the solutions, making this 
picture not free of certain assumptions. 

Assuming periodic boundary conditions (BCs) for N 
polarized electrons in a square of size D, one can ignore 
the constant term in He, the electronic density = N/D"^ 
and a — \/ ^/wnl. 



3 Square lattice model 

We now define a square lattice model of spacing s, size 
L = D/ s, nearest neighbor hopping element 



where = a/ as- The Bohr radius = h^/me^ is a 
measure of the GS radius of the hydrogen atom while the 
rydberg (IRy — /{2aB)) is its GS binding energy. Eq.|31 
tells us that the physics of a system of interacting electrons 
in the continuum does not depend on many independent 
parameters {h, e, m, the electronic density n^...) but only 
on a single dimensionless scaling ratio = a/os when 
TV — !■ cxD, aB characterizing the scale for the quantum 
effects. 

For the GS, if many electrons are inside the quantum 
volume a^, one gets the weak coupling limit (small rs) 
where one has a Fermi liquid. Though our theory could 
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and interaction strength 



2ms^ 



U = 



The lattice Hamiltonian Hi reads: 
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Fig. 1. Weak coupling limit: GS occupation numbers of a non 
interacting system in the reciprocal lattice. The Fermi-surfaces 
are sketched for increasing numbers of particles in a 36 x 36 
square lattice. At low fillings, the Fermi surface is almost a 
circle, while it becomes deformed at larger fillings. 



The operators cj (q) create (annihilate) a polarized elec- 
tron (spinless fermion) at the site j and (j, j') means that 
the sum is restricted to nearest neighbors, djy is the dis- 
tance between the sites j and j' in miit of s. 

The Hamiltonian (jS)) can also be written using the op- 
erators (dfc) creating (annihilating) a polarized electron 
in a plane wave state of momentum k: 

Hi = 4:Nt - 2t ^ (cos kx + cos ky) dl^dy^ 



where 



k,k',q 



. 1 cosqj 



(9) 



(10) 



The states of different total momenta K are decoupled. 
In the lattice units, IRy — U'^/At, gb = 2st/U and the 
ratio r, becomes: 



UL 



(11) 



4 Lattice effects in the thermodynamic limit 

4.1 Lattice threshold in the weal< coupling limit 

In the limit — > 0, the GS energy is mainly kinetic. This 
is a consequence of the Pauli exclusion principle, which 
pushes the Fermi energy of the non interacting system 



Fig. 2. Strong coupling limit: GS occupation numbers of a 
Wigner solid in real space. The hexagonal Wigner lattice of 
electrons becomes commensurate with the Lx L square lattice 
if one takes N — L = 30. 



to much higher values than the classical Coulomb energy. 
The GS kinetic energies of Hi and He become different 
when the Fermi surface is deformed by the lattice. Such a 
deformation can give rise to nesting effects with important 
consequences discussed in the literature 7,8,9. The Fermi 
wave vectors kp are given by At — 2i(cos k^ + cos ky) = cp 
for a square lattice of Fermi level eF^ instead of t{k^ + 
ky)'^ = ep for the continuum limit. 

Expanding the kinetic part of the lattice Hamiltonian 
for small wave vectors, one gets: 



k<kf 



(12) 



For the first term, one recovers the kinetic term of the 
continuum expansion Q): 



2~ / ^ ^ k^d^k 



k<kf 



(27r)2 2m 



N^h^TT 2 



{NRy), 



(13) 



while the lattice correction reads: 



Jj^ikt + kl)d'k (14) 



24ms2 

which becomes, using 

h^N^n^ 2 ttN 

AEi = ^ - -NRy- — . (15) 
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The condition for having the lattice correction AEi (Eq. 
115(1 smaller than the continuum kinetic energy (Eq. I13|l 
yields 

> r: « — . (16) 



This estimate of r* is only valid when ~^ 0, since it 
neglects the effect of the interaction. When one turns on 
the interaction, transitions from states below the Fermi 
surface to states above it (see Hamiltonian Q and Fig. 
nj take place. This smears the Fermi surface, giving an 
uncertainty Akp to kp- This uncertainty is evaluated in 
Appendix El One gets: 



(17) 



Since coa{kF + Akp) w 1 — (fci? + Z\fci?)^/2 when kp + Akp 
is small (say kp + Ak < 7r/2), one finds that the previous 
estimate oc s /as for r* is increased by an interaction effect 
oc {s/aBf" for small Tg. 



4.2 Lattice threshold in the strong coupling limit 

In the strong coupling limit Tg ^ oo, the GS energy is 
mainly classical (Coulomb energy) with weak quantum 
corrections. The electron lattice minimizing the Coulomb 
energy can be different if the square lattice of the model 
is not commensurate with the hexagonal lattice that the 
electrons form in the continuum limit. This is one ob- 
vious source of lattice effects when — > cxd. One does 
not discuss it, restricting the study to values of N and 
L where the two lattices are commensurate and give the 
same Coulomb energy. In this case, the lattice can never- 
theless change the vibration modes of the electron system 
around its classical electrostatic limit, when 7i — > 0. Let us 
consider the leading quantum corrections to the classical 
energy. 

In the continimm model, the GS energy per particle in 
rydbergs is given by Eq. [S] The first quantum correction 
AE^{rs oo) to the Coulomb energy (oc r~^) is given by 
the zero point energy of the harmonic oscillations of the 
electrons around their lattice points: 



AE'^ 



oo) 



A. 

3/2 ■ 



(18) 



In the lattice model, the classical limit i — + is not 

1/2 

described by the expansion in powers of valid in the 
continuum limit (Eq.[SJ), but by a lattice perturbation the- 
ory where the small parameter is /U . Examples of this 
lattice expansion valid when /U —+ can be found in 
refs. |lUlllllll2lll3lll4j . Its dominant quantum correction 
AEQ{rs oo) to the classical Coulomb energy comes from 
the term ANt of the lattice Hamiltonian JH)) . Expressed in 
rydbergs {IRy = U'^/it) per particle, this gives 



AE' 



oo) 



(19) 





m* /mo 


£r 


a*B (A) 


s(A) 


s/ag 


(1) 


0.19 


12 


33,2 


2.35 


0.071 


(2) 


0.07 


13 


100,0 


4.0 


0.04 


(3) 


0.6 


13 


12 


4.0 


0.33 


(4) 


10 


10 


0.53 


3.8 


7.16 


(5) 


175 


20 


0.061 


2.85 


46.7 



Table 1. Typical physical parameters for two dimensional sys- 
tems of charges of increasing effective masses created in var- 
ious devices: (1) Si-Mosfet, (2) n-doped GaAs - GaAlAs het- 
erostructure, (3) p-doped GaAs - GaAlAs heterostructure, (4) 
cuprate oxides exhibiting high- Tc superconductivity and (5) 
layered sodium cobalt oxides Naa;Co02. 



The next quantum corrections of order t'^ /U become neg- 
ligible in the limit t oc h — > 0. 

Let us consider a low density of electrons in a very 
large lattice, where one has the same hexagonal lattice 
with the same harmonic vibrations than in the contin- 
uum. The quantum corrections to the Coulomb energy 
and Vs are then given by Eq. ^1 and Eq. ^] respectively. 
If one takes the classical limit ?i — > in such a system, 
AEqIts) will reach the lattice limit AEQ{rs) which can- 
not be exceeded. This corresponds to coupling strengths 
where the harmonic vibrations of the electron lattice be- 
come so small that the discrete nature of the available 
space |12II13) becomes relevant. The lattice threshold r* 
is then obtained from the condition AEQ{rs) « AEQ{rs) 
and the continuum approximation is only valid if: 



N 



45 



1 



1/3 



2/3 



(20) 



Using the lattice spacing to the Bohr radius ratio s/ub, 
one gets that the lattice GS can be described by a contin- 
uum theory in the thermodynamic limit if 



> r* ~ 0.55 



-) 

as J 



4/3 



(21) 



and exhibits lattice effects otherwise. 



4.3 Two dimensional systems of different 

If the effective mass of the carriers is m* in a medium 
of dielectric constant Sr-, one must replace as by the cor- 
responding effective Bohr radius ajj. Some typical values 
of m*/mo, Sr, a*g, s and s/a^ are given in Table Q for 
two dimensional systems of charges created in various sys- 
tems. In general, the lattice spacing s is always of a few 
angstroms, and Sr ~ 10. But the carriers can be light or 
heavy, as indicated in Tabled 

When they are light, as in Ga-As heterostructures or in 
Si-Mosfets, s is small compared to a|j, and r* is small. The 
immobile ions can be modeled by a continuum uniform 
jellium, unless one reaches the very large densities where 
the Fermi surface becomes deformed. These densities are 
out of reach in today's semi-conductor field effect devices. 
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Therefore, if one numerically studies charge crystallization 
for those systems using a lattice model, J7/t, L and N must 
be taken such that = {U L / 2t) / [Vt^) > r*. 

In the cuprate oxides exhibiting high- Tc superconduc- 
tivity, the effective mass of the carriers is heavier, though 
not large enough for having s/a*g above 23, the value 
for which 0.55{s/a%)'^/^ « 37. The deformation of the 
Fermi surface by the lattice yields nesting effects which 
can be responsible for singlet d-wave superconductivity 
[TK] in those oxides. The lattice filling is large, reaching 
the limit where one can take a Hubbard model near half- 
filling. When one goes away from half-filling by chemical 
doping, those systems can reach a more dilute limit where 
one gets a quantum fluid which can be described by a 
continuum theory and where the scaling parameter be- 
comes relevant. 

There are systems with much larger effective masses, 
as those described by heavy fermions theories where s /a^ 
can exceed 23. One example is given jl6j by layered 
Lithium or Sodium Cobalt oxides: Lia;Co02 or Na2,Co02, 
where the effective mass of carriers can reach « 200, a 
value more familiar to /-band heavy fermions than to 
d-band metals, as discussed by Roger and Shannon|T71. 
Those systems look particularly interesting for the sub- 
ject of this study, since one should observe by increasing 
the carrier concentration a continuum-lattice transition 
for the electron (or hole) crystal, before having the quan- 
tum melting of this crystal in the lattice solid regime, to 
eventually obtain at higher densities a quantum fluid with 
a deformed Fermi surface. 

In Fig. 121 the regime of validity of the continuum ap- 
proximation is given in the (s/a|j, a/a*^) plane, for typical 
Ga-As heterostructures, cuprate oxides or layered sodium 
cobalt oxides Na2:Co02. The dashed line ly = N / L"^ — 1 
corresponds to systems usually described by the half-filled 
Hubbard model, when the spins are included. Increasing 
the density, one goes towards a lattice regime where the 
continuum approximation breaks down, giving rise to a 
lattice liquid {s/a*g < 23) or to a lattice solid {s/a^ > 23). 



5 Scaling in a lattice model with a fixed 
number of particles 

Let us define the lattice parameter suitable when N is 
constant. 

ri = ^^rsi2V^). (22) 

So far, we have considered the realistic case where N is 
varied by an electrostatic gate or by chemical doping in 
a lattice where L, t and U, and hence r/ are given. One 
can also study the lattice-continuum crossover by varying 
the lattice parameter n in a system of N particles. If N 
remains constant in a lattice where one varies ri cx r^, 
one gets the continuum limit and its universal scaling for 
small values of ri while the lattice becomes relevant and 
the continuum scaling breaks down for large rj. 

This can be seen from Hamiltonian (|5J), where the com- 
ponents of the two dimensional vectors k, k' and q can 



0.55 (s/a*) 




Fig. 3. Sketch of the lattice and continuum regimes in the 
{s / a*g , a / a*g) plane. In the non-shaded part, the lattice GS can 
be described by a continuum theory. The shaded part below 
the dashed line v = N/L'^ = 1 (2 with spins) is forbidden in 
a lattice model (more than 1 (2 with spins) electron per site). 
The line — a/a*g ~ 37 gives the density under which Wigner 
crystallization is assumed to occur in the continuum. The 
thick line 0.55(s/a|j)''/^ gives the lattice threshold r* above 
Tg ~ 37, while the dotted line gives oc s/a% + 0(s/a%)'^ 
below Ts ~ 37. The three double arrows correspond to typi- 
cal Ga-As heterostructures, cuprate oxides and layered sodium 
cobalt oxides Naa:Co02. Increasing the density, one goes from 
a continuum liquid or solid towards lattice regimes. 

take the values 0,2tt/ L, . . . ,2tt{L — 1)/ L. If the Fermi en- 
ergy is sufficiently small for having 4 — 2 (cos fc^: +cos ky) « 
for all the states below the Fermi surface, the kinetic 
energy reads 

ANt-2t (cosfc„ -|-cosfc,j,) w ^^p2, (23) 

k<kF i=l 

where k^ = (27r/L)pi with Px,Py << L for all k^ < kj?. 
Expressed in rydbergs {IRy — f/^/4t), the diagonal ma- 
trix elements of Hamiltonian Q due to the kinetic energy 
depend only on the lattice parameter rf. 

-^Ep' = ^ Ep'P?^- (24) 

i=l ' \i=l / 

The diagonal Coulomb matrix elements are given by 
N{N - l)C/y(q = 0) - UV{k, - k,,), while the off- 

diagonal terms cx U{V{'ki^ — k^, ) — V{ki^ — k'-, )), where 

V{q) is given by Eq.^| If U is small, only the off-diagonal 
terms with small momentum transfers (q = (27r/L)p with 
Px,Py << L) play a role. Transforming the discrete sum to 
a continuum integral, one gets for the interaction matrix 
elements: 

2L ^ di^o/L L L 
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2L 



Jo 



1 g27r2pr 



d(r,0) 



d'r=^I{p). (25) 



The integral /(p) is independent |12| of L when Px,Py << 
L and the Coulomb matrix elements of Hamiltonian ^ 
become also a function of the lattice parameter r; only: 



Ry. 



(26) 



For N fixed, assuming that N is small enough for 
avoiding deformed Fermi surfaces without interaction, the 
low energy levels depend only on the lattice parameter r/ 
when ri is small. The question is to determine the interac- 
tion threshold r^* above which the off-diagonal interaction 
terms begin to delocalize this GS to states of higher mo- 
menta, where 4 — 2 (cos + cos ky) ^ . When r; exceeds 
this r*, the lattice GS ceases to be a function of r; oc Ts 
as in the continuum limit. 



6 Lattice effects for a few correlated particles 



We propose three criteria giving similar lattice thresh- 
olds for an interacting system of N polarized electrons, 
which will be carefully studied when = 3 in the next 
section. The first one is a delocalization criterion in the 
Fock basis of lattice site orbitals. The second one uses the 
invariance of the persistent current when one varies the in- 
teraction strength in the continuum, an invariance which 
can be broken by the lattice. The third one is based on the 
limit imposed by the lattice to the zero point energy of the 
harmonic vibrations of an N electron solid, as previously 
discussed in the thermodynamic limit. 



6.1 Criterion 1: Delocalization in the Fock basis of 
lattice site orbitals 

Let us consider the system of N particles in real space 
instead of reciprocal space, in the limit t — where the N 
electrons are localized on N sites (see Fig. O and form 

states I J) = c]i---4nI'^^ °^ energy Ecoui{J)- As one 
turns on t, one can expect that the lattice becomes ir- 
relevant as each electron ceases to be localized on a single 
site. In analogy with the problem of a single particle in 
a disordered lattice, one can use the criterion first pro- 
posed by Anderson ^S]: delocalization takes place when 
the hopping term t between directly coupled sites becomes 
of the order of their energy spacing AE. This criterion 
was extended to interacting systems in many different 
contexts: onset of quantum chaos in many body spectra 
[Tnill9ll2Uj and in the quantum computer core quasi- 
particle lifetime and delocalization in Fock space |22II23| . 
In our case, the states become delocalized in the many 
body basis built from the states | J) when the matrix ele- 
ment {J'\Hkin\J) of the one body perturbation Hkin oc t 
coupling a state | J) to the "first generation" of states \ J') 
directly coupled to it by Hkin exceeds their energy spacing 
AEcoui = Ecoui{J') - Ecoui{J)- This gives t > AEcoui- 



Applying this criterion to the GS, one obtains r^* from the 
condition 

t^AEcoui, (27) 

where AEcoui is the increase of Coulomb energy yielded 
by the hop of one particle localized on the GS configura- 
tion to a nearest neighbor site when t = 0. When t exceeds 
AEcouii the GS is delocahzed on the J-basis, and hence 
on the lattice, and the lattice GS behaves as the contin- 
uum GS. 



6.2 Criterion 2: Persistent currents 

Since a continuum model is invariant under translations, 
the motion of the center of mass can be decoupled from 
the relative motions. Thus the continuum Hamiltonian He 
(Eq. jS)) can be decomposed in two parts, one related to 
the center of mass motion which is independent of the in- 
teraction, while the second one contains only the relative 
motions and hence the interaction. This has a very impor- 
tant consequence for the persistent current / driven by an 
enclosed Aharonov-Bohm flux in a continuum model: / 
is independent of ri and keeps its non interacting value. 
For having the topology of a 2d torus enclosing along 
the a;-direction, one takes the corresponding curled BC in 
this direction, keeping periodic BC in the y-direction. For 
a sufficient ri, the electrons form a Wigner solid and the 
small relative motions cannot feel the BCs. In this limit, / 
is just given by the center of mass motion, which is inde- 
pendent of ri , and hence coincides with its non-interacting 
value. This point remains correct for small ri, as it was 
proven for Id-rings J24"25','26' and observed for d = 2 
|24j . In contrast, since the previous decomposition into 
two parts does not necessary hold for Hi, I ^ I{ri = 0) 
for a lattice when Hi and He have different GSs. The de- 
cay of / above r* (small t/U at fixed N and L) can be 
evaluated jlllll2lll4| by the leading contribution (of order 

TV) if^ cx t{t/U)^-^ of the t/U lattice expansion. The 
value of ri for which 



i{u = 0) « 



(JV) 



(28) 



gives the Criterion 2 for (see Fig. (21 right). Instead of 
I{'P), one can prefer to use the Kohn curvature Ck = 
9^£'o/9<?^ evaluated at ^ = or the GS energy change 
AEo = Eo{<P = 0)- Eoi-P = <Po/2) where <?o is the flux 
quantum. To apply 'Po/2 corresponds to have anti-periodic 
BC in the x-direction. 



6.3 Criterion 3: Lattice limit for the zero point motion 
of an electron solid 

This is the criterion that we have already used in the qual- 
itative discussion of the thermodynamic limit, and which 
we consider when N is finite. When t/U 0, the leading 
correction to the Coulomb energy of Hi is 4Nt. Since the 
correction Eyib{ri) to the Coulomb energy coming from 
the zero point vibrational motion of the continuum solid 
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cannot exceed this lattice limit ANt, rf can be obtained 
from the condition 



(29) 



assuming that the values of the lattice parameters can 
yield a Wigner solid for ri < . 



7 Numerical study of three polarized 
electrons 

When one takes periodic BCs, a convention has to be cho- 
sen for the distance r (Hamiltonian 0) or djy (Hamil- 
tonian JHJ). For a finite square with periodic BCs, one 
possible definition is given by: 




"jj' 



where dx 



mm{\dx,L - \dx\)'^ + mm{\dy\,L - \dy\y, 

(30) 

Hereafter, we refer to 



Jx-Jx and dy =Jy-Jy 
the corresponding l/|(ijj'| repulsion as the periodic singu- 
lar Coulomb (PSC) repulsion, since it has a cusp when the 
interparticle distance djy has one of its coordinates equal 
to L/2. This cusp being unphysical, we introduce the pe- 
riodic regularized Coulomb (PRC) repulsion, defined from 



— \ I sm' 

IT 



L 



(31) 



which locally coincides with the PSC repulsion, but re- 
mains analytic for all values of d^y when s — > 0. The PRC 
repulsion is essentially equivalent to the Ewald repulsion 
obtained from the periodic repetition of the considered 
system. 

Defining djj' with Eq. 1311 we calculate the quantities 
used in the different criteria for iV = 3 polarized electrons 
on a square lattice. We give in Appendix [CI the same anal- 
ysis defining d^y with Eq.l^Olinstead of Eq.|^ The choice 
of the PRC or PSC repulsions, or of the repulsion obtained 
after Ewald summation is arbitrary for three particles in a 
toroidal geometry. Nevertheless, it allows us to check if our 
three proposed criteria for r* give consistent results when 
the long range form of the Coulomb repulsion is changed. 
The presented results extend to larger L previous studies 
of the case = 2 and TV = 3 given in Ref.[n| and Ref-fH] 
respectively. 

For t = 0, the configuration of particles minimizing 
the PRC Coulomb energies is given in the inset of Fig. 
01 The values oi L = 6,9, 12, 15, 18, . . . yield a diagonal 
Wigner molecule shown in the inset which is commen- 
surate with the square lattice. Moving one particle by a 
single hop increases the Coulomb energy Eq — {^/GUn)/ L 
by an amount 



Coul 



7V2tt^U 
12V3i3 



(32) 



when L is sufficiently large. 

For [7 = 0, the GS energy is given by ^0(0) = 12<-8i- 
4tcos(27r/L) for periodic BCs and becomes Eo{'Pq/2) = 



Fig. 4. Energy ratio Fn=3{L, U,t) as a function of n — UL/t 
given by the PRC repulsion for L = 6 (>), 9 (□), 12 (o), 15 (A), 
18 (o). The dotted-dashed line gives the behavior 0.2327y^ 
(harmonic vibrations of the continuum Wigner molecule) and 
intersects the limiting dashed lines 12t/{At ~ it cos{2n / L)) at 
the r*{L) corresponding to criterion 3. Inset: A GS configura- 
tion when t = and L = 24. 



12i — St cos(7r/L) — 4i cos(37r/L) when one twists the BC in 
the x-direction. The difference AEq = Eo{^o/2)-Eo{0) ~ 
— Mn'^t/L'^ when L ^ 00. When t/U is small, AEq can 
be calculated at the leading order of a ^/{/-expansion 
for A^ = 3. This gives when L is large: 



hm AEo « — ^ ; hm AEq « -— 33 

n-.o L2 n^oo L^AEcoui 



where AEcoui is given by the Eq. Using these expres- 
sions, one obtains from the two first criteria: 



(34) 



where a = 4 for the PRC repulsion, the constant A slightly 
depending on the taken criterion. 



7.1 Persistent currents 

We now present numerical results obtained using the 
Lanczos algorithm, using Hamiltonian (|5J) and considering 
the sub-space of total momentum K = 12 for periodic 
BCs (no applied flux). 

Our system has the topology of a 2c? torus. To en- 
close an Aharonov-Bohm flux along the x-direction, one 
takes the corresponding curled BC in this direction while 
the BC in the y-direction remains periodic. To apply half 
a fiux quantum (^ = ^o/2) is equivalent to take anti- 
periodic BC along the x-direction. In Fig[Sl the increase 
-E'o(^) — Eo{0) of the GS energy Eq is given as a func- 
tion of 'P/'pQ for different values of r; using a 18 x 18 
square lattice. When r; is small, the curves coincide. This 
is the continuum regime where the persistent current is 
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Fig. 5. GS energy Eo{<P) — Eq{<1> = 0) as a function of the 
enclosed dimensionless magnetic flux ^/^o for A'^ = 3, L = 18, 
PRC interaction at r; = 6 (o), 60 (□), 600 (o) 6000 (A) and 
60000 (<). 



AE„(r,)/AE„(r=0) 



0.1 - 




Fig. 6. Dimensionless change AEQ{ri) / AEi){ri = 0) of the GS 
energy when the longitudinal BC is twisted for L = G (□), 9 
(o) ,12 (A), 15 (<), 18 (o) and A' = 3 as a function of n (PRC 
repulsion) . 



independent of the interaction. When r/ is large, the in- 
crease Eo{(I>) — Eo{0) becomes weaker. This is the lat- 
tice regime where the persistent current decays as the 
interaction increases. One gives in Fig|H| the dimension- 
less change AEo{ri) / AEo{ri = 0) of the GS energy when 
the BC is twisted in the a;-direction for increasing val- 
ues of L. One can see the two limits given by Eq. 123 
AEo{ri)/AEo{ri = 0) w 1 in the continuum limit, fol- 
lowed by a decay when r; exceeds the lattice threshold 



B = {V6e'^TT)/{2AD^m). The zero point vibrational energy 
is then given by: 

Eyib{rs,N = 3) = h{LUi + LOt) 

in rydbergs where /3 = 3/2, with Vs — n/(2V37r) for iV = 
3. 



7.2 Harmonic vibrations of the continuum Wigner 
molecule 

For the third criterion, one needs the zero point vibra- 
tional energy of the Wigner molecule that the three par- 
ticles form when r/ is large, but smaller than r^*. This can 
be calculated in the continuum limit, using for N — 3 the 

1/2 

same expansion in powers of Ts than those used in Eq. 
1^1 for N oo. We summarize the main points, the de- 
tails being given in Appendix ^ In the continuum, the 
Hamiltonian He can written as the sum of two decoupled 
terms. Denoting R = (^^ ri)/3 the coordinate of the cen- 
ter of mass, the first term reads Hcm — (^^/6'ti)V^ and 
corresponds to the rigid translation of the molecule while 
the other term contains the relative motions and the in- 
teraction. For a Wigner molecule, the second part can be 
simplified and expressed in terms of the normal coordi- 
nates suitable for describing the small vibrations around 
equilibrium. 

The PRC repulsion is harmonic around equilibrium, 
and the three particles form a diagonal chain as indicated 
in the inset of Fig. ^ when L/3 is integer. One gets four 
decoupled harmonic oscillators, two corre sponding to a 
longitudinal mode of frequency uji = the two others 

being a transverse mode of frequency ujt — V8i?, where 



7.3 Scaling of the ground state energy 

From the GS energy Eo{L, U, t) of K = 0, and for a given 
value of N, we define the dimensionless ratio Fm[L, U,t) 
by: 



FN{L,U,t) 



Eo{L,U,t)~ EoiL,U,t^O) 
Eo{L,U = 0,t) 



(36) 



This ratio gives the change of the GS energy from the 
Coulomb energy due to the quantum effects, divided by 
the GS energy without interaction. 

The results for the PRC repulsion are shown in Fig. 
H For i = 0, the values of L = 6,9,12,15,18 are 
commensurate with the period of the diagonal Wigner 
molecule shown in the inset. This gives the same classical 
Coulomb energy for the lattice and the continuum when 
t ^ 0, eliminating a trivial source of lattice effects. When 
FN^3{L,U,t) is plotted as a function of r;, the different 
functions F/v=3(L, U, t) scale without an observable lattice 
effect up to the r*[L) exactly given by Criterion 3. Using 
Eq{L,U = 0,t) = 12t-8t-4tcos{27T/L) one can see that 
the numerical results coincide with the analytical result 
i^7V=3 = 0.2327-^77 implied by Eq.|211for intermediate val- 
ues of r; where one has a continuum Wigner molecule. The 
function FN=3{L,U,t) saturates to 4:Nt/Eo{L,U = 0,<) 
above rKL), as indicated by the dashed lines. 
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Fig. 7. Energy ratios Fn{L, U,t) using tlie PRC repulsion for 
iV = 3 (i = 18 solid line) and = 4 with L = 6 (o), 8 (□), 10 
(A) as a function of — ri / (2\/ A'"7r). (PRC repulsion). 

8 Effect on the scaling function when 
varies 

In Fig. [71 a small change of the scaling curve F]y{rs — 
ri/ {2^/ N-k) can be seen when a fourth electron is added, 
accompanied by the expected breakdown of the scaling 
behavior above the corresponding r,* for = 4. When 
N oo, should converge towards a thermodynamic 
limit depending only on . Unfortunately, a study of this 
convergence is out of reach of a numerical approach using 
exact diagonalization. 



9 Conclusion 

We have studied the lattice effects upon an interacting 
system of polarized electrons in two dimensions. We have 
first considered the case where the number N of polarized 
electrons is increased in a square lattice of large size L and 
of fixed parameters U = e^/{ers) and t = / {2m* s^) . 
This corresponds to semi-conductor field effect devices or 
layered oxides where the number of carriers can be varied 
by an electrostatic gate or by chemical doping. Starting 
from an empty lattice, one has a continuum regime and 
its universal scaling if one uses the parameter Tg, until 
the carrier density n* = l/(-y7rr*a|j)^ is reached. At this 
density, the continuum approximation with its universal 
scaling breaks down, as sketched in the phase diagram 
given in Fig. O This lattice threshold r* takes place in a 
quantum fluid phase if the carriers are light and the ratio 
s/a*g < 23. We have studied more particularly in the re- 
maining part of our manuscript the case of heavy carriers 
where s/a*^ > 23, for which the continuum approximation 
breaks down in the crystalline phase. We have pointed 
out the limit imposed by the lattice of positive ions upon 
the zero point motion of the electron lattice. We have ne- 
glected the obvious problem coming from the incommen- 
surability of the two lattices, and focus our attention to 



the commensurate case. The studied lattice effects are in- 
dependent of this commensurability issue, which should 
matter at large lattice fillings. 

In the second part of this manuscript, we have studied 
the role of a lattice for a fixed number = 3, 4 of po- 
larized electrons. The lattice-continuum crossover is then 
obtained by varying the lattice parameter r; = (UL)/t. 
ri and Vs play the same role when N is fixed, since 
Ts — ri/{2\/nN). The continuum approximation is valid 
and there is a universal scaling when one uses the pa- 
rameter ri as far as ri does not exceed a lattice threshold 
T;* which has been determined from three criteria. One 
of them was based on the behavior of the persistent cur- 
rent I{ri) driven by an enclosed Aharonov-Bohm flux. For 
n < f*{L), I{ri) = I{ri = 0) while /(n) decays above r* 
and ri ceases to be a scaling parameter. For a finite number 
of particles, one goes from the continuum regime towards 
the lattice regime through a smooth crossover. The ques- 
tion to know if this smooth crossover does not become 
sharper when iV ^ oo, to give rise to a true quantum 
transition is an interesting issue which we postpone to a 
following study. 

Z. A. Nemeth acknowledges the financial support provided 
through the European Community's Human Potential Pro- 
gramme under contract HPRN-CT-2000-00144 and the Hun- 
garian Science Foundation OTKA T034832. 



A Weak interaction correction to r* 

When Ts is small, the main effect of the interaction is 
to smear the Fermi surface, giving an uncertainty Akp to 
kp, such that one expects to have the continuum behavior 
when kp -\- Akp, and not only kp, is small. To evaluate 
Akp, we assume that the low excited states only become 
occupied at low Vs- The first excitation energy reads: 



N 

AE,,o=ty(kl 



i=l 



= 2i(k'2 



kp) w AtkpAk 



(37) 



where the factor 2 comes from momentum conservation 
and k' is the wave vector of an empty state above the non 
interacting Fermi surface. This gives us the relation: 



AEf 
'Ak'p 



Atkf. 



(38) 



The Fermi energy uncertainty AEp can be estimated from 
the spreading _r of a non interacting level, when one turns 
on the interaction. Using Fermi's golden rule, one gets 



F : 



27r|i?o,i" 



n(kp) 



(39) 



where the matrix element of interaction coupling the GS 
to the first excited state |i7o,i| = 2C/(U(q) - V{2kp)) « 
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2C/y(q) reads 
U 



J 



V 

T 



cos 2lTX 



lU 
T' 



(40) 



/ is a constant equal to 1.029 for the PRC repulsion. The 
number of states n{kF) on the Fermi surface is equal to 
V4iV7r and AEiq = '^tkp\c{\ where |q| = {2tt / L) is the 
smallest momentum for an excitation. 

One eventually gets for the Fermi energy uncertainty 



r : 



it 



which gives for the Fermi momentum uncertainty 



UI 
T 



L 



32VttN 



The condition 



is satisfied if 



kp + Akp < 



TT 



16 



2 

4 s 
TT as 



(41) 

(42) 
(43) 
(44) 



When rs is not too large, the continuum theory is valid 
if Vs > r* with a threshold r* having a small correction 
cx (s/as)'^ driven by the interaction: 



4 s 

TT Qb 



4/2 



(45) 



The constants in the expression of r* depend on the used 
criterions for neglecting lattice effects (for instance A/tt 
comes from the condition H43|l ). 



B Zero point energy of a continuum Wigner 
molecule for = 3 

For three spinless fermions on a continuum square domain 
of size D with periodic BCs, the continuum PRC repulsion 
reads 



Vir) 



Dv/sin 



2 



D 



sm 



2 r-yTT 



(46) 



If D is large enough, the GS is a "Wigner molecule" of 
delocalized center of mass, but of quasi-localized inter- 
particle spacings for minimizing the Coulomb energy. For 
a certain center of mass, the molecule of lowest Coulomb 
energy with the repulsion l^Hjl consists in putting the par- 
ticle coordinates at ri = (0,0), r2 = {D/3, D/3) and 
r3 = (— D/3, ~D/3). This configuration has the Coulomb 
energy 

Ecoui = V6e\/D. (47) 

The particles forming this molecule vibrate around the 
equilibrium positions. This motion is an harmonic oscilla- 
tion if the amplitude of the vibration is small. To describe 



this harmonic motion, one expands the pair-potential H46|l 
around the equilibrium distance Tq — (-D/3, Z3/3) up to 
the second order: 



+ (•••) 




D 




-|)+0(r3), (48) 



where the missing term (...) is the first order contribution 
which disappears after summing over all the pair poten- 
tials. The expansion H48|l becomes: 



V{r)'. 
where 



3D 



(...) + (r-ro) 



^)(r-ro). (49) 



(50) 



72 D3 
and B = 3 A/7. 

The three particle Hamiltonian with the expanded re- 
pulsion becomes He ~ Ecoui+Hharm, where the harmonic 
part is: 



harii 



2m 



(V? 



V^) + XMX. 



(51) 



The vector X = {xi, yi, X2, 2/2, a^3, ys) is composed of the 
6 relative coordinates and the 6x6 matrix M is given by: 



M = 



( 2A 


2B 


-A 


-B 


-A 




2B 


2A 


-B 


-A 


-B 


-A 


-^A 


-B 


2A 


2B 


-A 


-B 


-B 


-A 


2B 


2A 


-B 


-A 


-A 


-B 


-A 


-B 


2A 


2B 


\-B 


-A 


-B 


-A 


2B 


2A ) 



(52) 



Diagonalizing Af, one obtains the normal modes of the 
harmonic oscillations while the eigenvalues of M give their 
frequencies. One obtains 

— Two eigenvectors of eigenvalue 0. 



Xi = -^(1,0,1,0,1,0)-X, 
X2 = -^(0,1,0, 1,0, 1)-X. 



(53) 



This zero frequency mode corresponds to the transla- 
tion of the center of mass of the molecule. 
— Two other eigenvectors of eigenvalue lO-B, correspond- 
ing to the longitudinal mode (vibration parallel to the 
axis of the molecule). The normal coordinates can be 
taken as: 



X3 - ^(i,i,-i,-i,o,o)-x, 

X4 - 4^(1, 1,1, 1,-2, -2) -X; 



(54) 
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Two eigenvectors of eigenvalue 4i?, corresponding to 
the transverse modes. The normal coordinates can be 
taken as: 



X5 = ^(i,-i,-i,i,o,o)-x, 

X6 = 4^(1,-1,1,-1,-2,2).X. 



(55) 



Using these normal coordinates, the Hamiltonian H51|l 
becomes a decoupled sum of two harmonic oscillators: 



harm = " E + ^^^ixl + xl) + + xl) . 



(56) 

For a GS of total momentum K = 0, there is no motion of 
the center of mass, the GS wave-function does not depend 
on Xl ^-iid X2 and can be factorized as: 



,X&) = '^aL{xz)VoL{Xi)VoT{Xb)VoT{x&) (57) 



where L, T refers to the transverse and longitudinal modes 
and tfQ to the ground state of an harmonic oscillator: 



<po(a;) 



1 



of length = (h^ / (rn^uj^)) . One eventually obtains for 



exp - 



(58) 



the GS energy with the expanded pair potentials: 

Eq - Ecoul = Hi^T + LdL)\ 
(^L = 

UJT = 

and using the expression of B: 

Eq — Ecoul 




(59) 
(60) 

(61) 



20VQh'^e^Tr3 / 8\/6 ;i^e%3 



24 £137 



24 D^m 



\/2)i 



(62) 



For the energy ratio Fn=3{L, U, t), using for the kinetic 
energy in the continuum limit Eo{L, U = 0,t) = Sn'^t/L'^, 
one gets the behavior numerically obtained from the lat- 
tice Hamiltonian Hi and shown in Fig. ^ for intermediate 



Fairs) 



V5 + V2 VGtt^ fUL 



87r2 



3 V t 

1/4 



= 0.5764^^ = 0.2327Vn- 



(63) 




Fig. 8. Energy ratio Fn=3{L, U,t) as a function of n — UL/t 
using the PSC repulsion for L = 6 (>), 8 (□), 10 (o), 12 (A), 14 
(o), 16 (-I-), 18 (<). The dotted-dashed line gives the rf^^ be- 
havior due to the vibrations of the continuum Wigner molecule. 
Inset: a GS configuration when t = and L = 24. 



C Lattice threshold r/ 
repulsion 



using the PSC 



With the distance djj/ defined by Eq. we have previ- 
ously studied the validity of a continuum approximation 
for a lattice model of three polarized electrons. In this ap- 
pendix, we revisit the same issue defining djj' from Eq. 
1^ instead of Eq. El Let us calculate the quantities used 
for the three criteria when one uses the PSC repulsion. 
For t — the Wigner "molecule" minimizing the PSC 
Coulomb energy has the triangular shape shown in the 
inset of Fig. |S1 instead of the linear shape shown in the 
inset of Fig.^ Moving one particle by a single hop in this 
triangular molecule increases the PSC Coulomb energy by 
an amount 

(64) 



^Ecoul 



L2 



,(Pi?,C) 



when L is sufficiently large, instead of the AE)j^^i 
U/L^ given by Eq. (23 

For the energy change AEq, one obtains the same ex- 
pressions as in Eq. but with AEcoui given by the Eq. 
|M] instead of Eq. jJH The two first criteria gives 



^(i) = AL", 



(65) 



where a = 3 for the PSC repulsion, instead of a = 4 for 
the PRC repulsion. 

When one takes the PSC repulsion, the three relative 
distances at equilibrium are precisely r — (L/2, L/2), r = 
(0, L/2) and r = (L/2, 0) respectively when L is even. The 
potentials v{5v) felt by the electrons around their equilib- 
rium positions are singular and instead of the analytical 
expansion of w ((5r ) , one has v{5t) « Ci | Sr^ | + C2 1 5ry \ , 
where Ci and C2 depend on the equilibrium positions 
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and are cx / — Ujl?. For a single particle in a Id- 
potential vix) = C|a;|, the GS energy e can be approxi- 
mated by t/B^ -I- CB where B is the GS extension and 
is given by de/dB = 0. This yields B (x (C/t)^/^ and 
e DC {U'^t/L^y^^ . Since the 2d-potential w((5r) is separa- 
ble, one eventually finds: 



remain valid independently of the used definition of the 
Coulomb repulsion in the periodic square lattice, 



(66) 



in rydbergs where /3 = 4/3. As one can see, the PSC re- 
pulsion gives a higher exponent /3 when N = 3, which 
is inconsistent with the usual expansion [2] in powers of 

— 1/2 

Ts first proposed by Wigner. 

Using Eg. 1661 one gets from Criterion 3 given by Eq. 
IMI again, but with a = 3 for the PSC repulsion instead 
of a = 4 for the PRC repulsion. The PSC repulsion is 
somewhat unphysical and leads to stronger lattice effects, 
but provides an interesting check of the validity of our 
theory: The changes of AE^^^f'' and e''^^^^ are such 
that the different criterions give thresholds which are 
consistent. 

The dimensionless energy ratio i^jv (£, C/, t) for the PSC 
repulsion is shown in Fig. |H1 for even values of L, where 
the GS is a triangular "molecule" shown in the inset when 
t/U 0. Again the curves scale up to the onset r;*(L) 
given by Criterion 3. But the PSC repulsion gives rise to 
a different onset r^i^L) than the PRC repulsion for = 3, 

2/3 

since at intermediate ri one has Fm=3 csc for the PSC 

1/2 

repulsion, and not oc as for the PRC repulsion. 

Does this difference remain for larger values of A^? In- 
deed the contribution of pairs ij having the coordinates of 
their spacings dij close to D/2, and responsible for the r1^^ 
behavior when dij is defined by Eq. 1301 becomes a surface 
effect oc A^ compared to the bulk contribution c>c A^^ of the 
remaining pairs, yielding AE^^^i w AN jl? -|- BN"^ jl? , 
where A and B are constant. For a fixed L and increas- 
ing Af, AE^^^i BN'^/L^ and following Criterion 1, the 

1/2 

conventional expansion for should be valid for the 
PSC repulsion too. Therefore, large periodic square lat- 
tices should exhibit a behavior independent of the choice 
of the long range part of the Coulomb repulsion when A^ 
becomes large. Another possible choice is the Ewald repul- 
sion obtained after summing over all the electrons present 
in the infinite repetition of the same finite square in the 
X and y directions. For a small number A^ of electrons 
in a periodic square, these definitions are somewhat ar- 
bitrary. But to reach the thermodynamic limit, the PSC 
repulsion is less appropriate than the PRC or Ewald re- 
pulsions, since it gives larger finite A^ effects. 

Nevertheless, the following relations for the lattice 
threshold, the continuum zero point energy of the crys- 
talline oscillation and the characteristic scale of the 
Coulomb energy respectively: 



ri (X L 
Evib oc r^'^ 
AEcoui cx UL-^ 



(67) 
(68) 
(69) 



a = 7 + 1 ,a = 



2-/3 



(70) 



between the exponents. 
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